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Abstract 

Simplified solutions of the Cox- Thompson inverse quantum scattering method at fixed energy are 
derived if a finite number of partial waves with only even or odd angular momenta are contributing 
to the scattering process. Based on new formulas various approximate methods are introduced which 
prove applicable also to the generic scattering events. 
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1 Introduction 

The Cox- Thompson (CT) inversion method [I] has a number of useful properties. As input it requires a 
finite set consisting of N phase shifts. As output it produces potentials which possess non- vanishing first 
momentum, J rV(r)dr ^ 0, and finite values at the origin, |V(0)| < oo [2]. These properties make the CT 
. method an attractive constructive procedure which can be applied to experimental phase shifts in order 

to determine interaction potentials between composite quantum systems. However, the CT method leads 
to the solution of equations which are of nonlinear character, involving the inverse of a matrix containing 
the unknown quantities themselves. The solution of the CT equations is therefore numerically difficult 
to perform and, sometimes, only successful with deployment of sophisticated nonlinear solvers such as 
the simulated annealing method [5]. 

In this paper we investigate the possibility of making the CT equations simpler and introduce some 
alternative forms which follow from the analytical structure of the method. We shall derive equations 
for calculating the unknown generalized angular momenta, L's, from the finite input set of phase shifts, 
Si's. The characteristic property of these equations is that they do not involve any complicated matrix 
inversion. The simplification is only possible if the even angular momenta can separately be handled 
from the odd ones. Physically such a situation occurs in certain case of identical bosonic (fermionic) 
collision when only even (odd) partial waves are contributing to the differential cross section and, thus, 
only phase shifts belonging to even (odd) angular momenta can be derived from experiments. We note 
that according to references [3j 0] , the even (odd) partial wave phase shifts alone would result in the 
same (unique) potential if an infinite subset of them would be used as input and a proper constructive 
procedure would be known. In case of finite set of phase shifts, however, when also even and odd partial 
waves contribute generally to the scattering, we can still use our simplified forms as approximations, and 
we shall provide various approximate examples for such usages. Finally, we also investigate the possibility 
of using the so called one-term solution (or total decoupling case) approximation L = I — 26 g/ 'it, which 
can be obtained from the CT equations with N = 1. 

In the next section we collect all the equations necessary to outline the CT method. The simplified 
equations valid in case of even or odd angular momenta are derived in section 3 and the examples for the 
various applications with and without approximations are presented in section 4. Section 5 is devoted to 
the conclusion. 
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2 Cox-Thompson method 

Let us introduce first the dimensionless quantities, x = kr for the distance, and q(x) — YlSiLMl f or the 
potential where E = h 2 k 2 /2m denotes the scattering energy, m the reduced mass, and k the wave number. 
The inverse potential corresponding to the radial scattering Schrodinger equation reads as follows 

2d K(x,x) 

q{x) = — x Tx—r- (1) 

where the quantity K{x, x) is the diagonal part of the transformation kernel defined by a Gel'fand- 
Levitan-Marchenko (GLM) type integral equation 

K(x,y) = g(x,y) - [ dtr 2 K(x,t)g(t,y), x>y. (2) 
Jo 

In the CT method we take the following separable expansion for the input kernel: 

g(x,y) = ^2nM x <) v i{ x >), \ ^ \ = \ ™ \( x ,v) ( 3 ) 
les l > ) l J 

where 7; are expansion coefficients (not required to be determined) , uj (vi ) denote the regular (irregular) 
solutions of the free radial Schrodinger equation and the summation set S consists of the physical angular 
momenta {/}. 

To solve the GLM equation Cox and Thompson [T] introduce another separable ansatz for the 
transformation kernel, in the form of a sum over an artificial angular momentum space L G T 

K(x,y)=Y,M x W{y) (4) 
let 

where A^s are unknown expansion functions and T is interpreted as a set of unknown shifted angular 
momenta to be determined under the constraint that it has the same number N of different elements as 
S does, i.e., AT = \S\ = \T\ and S n T = 0. 

By inserting equations ([3} and (gj) into the GLM equation @ and using the linear independence of 
the regular and irregular free solutions, one obtains equations for the determination of the expansion 
functions Al(x) as follows 

Here, the only unknowns are the set T with elements L, and W denotes the Wronskian defined by 
W[a, b] = ab' - a'b. 

In order to determine the set T one makes use of the Povzner-Levitan representation of the radial 
scattering wave function which reads as 

ifo(x) = ui(x) - I dtt- 2 K(x,t)ui(t), I G S (6) 
Jo 

whose asymptotic form containing the input phase shift data, takes the form 

Here, Bi is a normalization constant and we have defined the asymptotic expansion functions A\[x) = 
Al(x — > 00) which can be calculated from the asymptotic version of equation ([5]) which is given by 

E4W/f#T = -»(x-!^ ^G^. (8) 
iW Z(/ + l) -L(L + 1) y 2 h w 

Using the last two equations, ([7]) and ©, one can easily derive the following equations for the determi- 
nation of the unknown L's from the input phase shifts, Si's 

Sl = - or tan(^)= + — , l&S (9) 
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with 



where Si — e 2lSl and the "shifted" reactance matrix elements are defined as 

Kf= [MMM-^Ure^ 1 - 1 '^ 2 , leS, (10) 

LeTMes 

M cos j lL - L(L + 1) - 1(1 + 1) \ cos ((I - J ' Uij 

The appearance of the inverse of the matrix M cos containing the unknown L's is the problem that makes 
the solution of the nonlinear equations © especially hard. 

If one succeeds to solve either one or both of the highly nonlinear equations © and thus finds the set 
T, the CT inverse potential qcr(x) can easily be obtained by employing equations ©, (jl]), and ((!]). 



3 Simplified solutions and approximations 

In this section we present simplifications to equations © which can be used if only even (odd) partial 
waves are arising during the collision. Otherwise the simplified equations can be employed to construct 
different approximations which will be discussed in separate subsections. 



3.1 Even (odd) angular momentum treatment 

Let us differentiate equation ([5]) twice with respect to the variable x. Then we arrive at the following 
equation 

= -Al(x), (12) 

which has a periodic solution as 

A*L(x) = a L cos(x) +b L sin(x). (13) 

Now, by inserting this solution ([13")) into equation © and taking into account the independence of the 
sine and cosine functions, one gets the following two equations 

yf^l cos ((/-£)§) Jco S (Zf)\ 
^1 b L j L(L + l)-l(l + l) = { sin (if) f< lES - 



LET 

Consider the decomposition 



S = S e U S 



where S c and S contains, respectively, the even and odd elements of S. Instead of (14) we consider two 
systems: 

>p f a L \ cos(q-L)f) _ f cos (If) 1 

1 br [ t.(t.+ -\\-ui + 1 sin p fc ° c 

and 



L — L b L f L(L+1)- 1(1 + 1) \ sin (if) 



v J i cos(q-L)f) r cos (if i L 

2- 1 b L f L(L+l)-l(l + l) { Bin(Zf) /' <e *°' 



LeT 

where |T e | = |£ e |, |T Q | = |5 Q | and T e n 5 e = 0, T Q n 5 Q = 0. These systems have the solutions 

fli = rf / r / r rri — TTTTi ZTvi rT^v &i = °' X e Tc ' ( - 15 ^ 

lli'STeXCi}^^ + 1 ) ~ L \ L + 1 )J C0S ( L f ) 

a L = 0, b L = — 77m7T\ r7777~TTvi ■ tT*\ > l ^ t o, (16 

LLL>£T \{L}\ L ( L + 1) - L \ L + 1)) sm ( L f ) 

respectively. In case of T n T ^ the formulae (fT5|) and (fTB]) may assign different values to the same a l 
and 6l but this is not a real ambiguity because we always use separately the solution vectors (fT5)l and 
(fl6l) in our calculations. 



and 
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Now, by using the explicit expressions (fT5j) in equations (fT3|) and ([7]), one obtains the final solution 
to the CT method as 

( ° = " S. Orf + 1) - ^ + D) ( 2) < ^ (17) 

valid for the case of even Vs. Similarly, using equations (fT6|) we get the solution to the CT method as 

tan(5 ° = S rww^ + 1) - l>(l> + d) cot ( L 2) ' leS ° (18) 

which are valid in the case of odd Vs. These equations determine the unknown sets T c or T a of shifted 
angular momenta L and, thus, replace either of the corresponding generic solutions ©• 

Notice the simplified structure of the nonlinear equations ([IT]) and (|T8|). compared to equations ([9]). 
While equations © contain an explicit matrix inversion of a matrix involving the unknowns of shifted 
angular momenta, L's, formulas (1171) and (1181) do not require such a nonlinear operation. They 'only' 
contain products and the tangent (cotangent) operations and are thus presumably easier to be solved for 
the sets T e or T a , if the respective input phase shifts are given. 

Finding the sets T e or T Q , the corresponding potentials q e (x) or q Q (x) can be obtained similarly as in 
the general case, by employing equations ©, ((!]), and (TTJ). 

3.2 Equivalence of solutions (JSj) and flUD or (USD 

By explicit calculation one can check the equivalence of equations © and (|T7|) or (fT8"|) for the special cases 
of even or odd Vs. For the case of either even or odd Vs, the relation /C ; + = /C ; ~ holds. Now, specifying 
ourselves to the even I case only, I £ S e , the second of the general solution (J9j> can be written as 

tan(<y= E [M s m}lL[M-l} Ll '(-f- n/2 , left. (19) 

Using equations flHJl, and (|T3|). we get the expression 

a L cos(x) = ^ [M^]/* cos(x - L £ T c (20) 



i'es c 



which simplifies to 



o £ = S^lwH' 72 ' L G T c . (21) 

By multiplying both sides of this equation by (— )^ 2 [Af s i n ];i, and performing the sum over L's, one may 
write 

a L [M sin ] lL (-) l/2 = ]T [M 8in ] a [M-i]^(-) ( ^ 0/2 , l£S e . (22) 
LeT e LeT a ,i'es c 

According to equation (|T9")) . the right hand side is already equal to tan(<5j), and, by noting that the 
matrix M s ; n on the left hand side can be written, on account of (fTTj) . as [M s i n ]iL = — (— )'^ 2 sin(£?)/(i(i+ 
1) — Z(i + 1)), one gets the formula 

tan ^ = - E ^Z(LTTr|Ti)' (23) 

which is the same as equation (fT7]) if one takes into consideration the solution (|15D for the coefficient 
a,L,L £ T e . 

A similar procedure can be applied to proving equivalence of equations © and (fT8")) for odd Vs 
3.3 Approximations 

Equations (TTTJ and (TTS)) can be used to derive inverse potentials only in the case when either even or odd 
partial waves are arising during the collision process. In other words, only identical bosonic or fermionic 
scattering can be treated by the simplified equations (|T7l) and (|T8|). However, the simplified equations 
offer several possibilities to introduce various approximate treatments of the general scattering case. The 
type of approximations will be classified according to the level it is applied to. 
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3.3.1 Potential approximation, A 

If the general equations © can be solved by neither of the nonlinear solvers at hand, one may try to 
assess the inverse potential by solving the simplified equations (TIT)) and (1T51) for the sets T e and T Q . Then, 
separately constructing the corresponding potentials q e (x) and q (x), one simply adds them together to 
get an approximation of the interaction potential, qA(x) — Qc(x) + q {x). 

3.3.2 T-set approximation, T 

One may try to approximate the set of the shifted angular momenta themselves. By unifying the two 
sets obtained by solving equations (TIT)) and (fT5|) . one gets the T-set approximation T a = T U T D . Using 
this approximate set T a in conjunction with equations dSJ , d4j) , and (HJ, one gets the approximate inverse 
potential qr(x). 

3.3.3 One-term approximation, L 

If the collision is dominated overwhelmingly by a single partial wave (as in the case of resonance scattering) 
then the equations © are to be solved at N = 1, and this results in the simple expression L = I — 25i/ir 
for the shifted angular momentum, assuming that the Ith partial wave is dominating. If however this is 
not the case, one still may try to use the approximate expressions 

L a = l- 25,/tt (24) 

to form an approximate set Tl. Using this approximate set Tl in conjunction with equations ([5)1, Q, 
and (TT]), one gets the approximate inverse potential denoted by qi,(x). 

4 Examples 

In this section we first apply the simplified solution (IT71) to calculate an effective potential related to 
phase shifts derived from bosonic collisions. Then, two exploratory calculations follow demonstrating 
applicability of different approximations introduced in the preceding section. 

4.1 Effective 87 Rb + 87 Rb potential at E = 303 ^iK 

In this subsection we recalculate one result of [7] where effective Rb-Rb inter-atomic potentials have been 
derived from ultracold Bose-gas collision data. 

The inverse calculation is performed by using equation p7|) at E = 303 /iK. The corresponding phase 
shifts <5° 1Ig measured (in rad) are as follows: —1.287, 1.635, 0.005 related to partial waves with I = 0, 2, 4, 
respectively. 

The resulted inverse potential Vqt (r) is plotted in figure [T] To control the procedure we recalculated 
the phase shifts from the inverse potential and the difference A CT = |<5° ns — <5; CT | between the phase 
shifts are as follows: 0.014, 0.042, 0.004 for I = 0,2,4, respectively. 

The result shown in figure [T] is identical with that obtained by [7j using the generic CT procedure. 
Indeed we also calculated the inverse potential by using equation ©, and obtained the same result as 
depicted in figure [1] In other words, the equivalence of equation (|TT)) and (JSJ) in case of even l's has been 
numerically verified too. 

4.2 Gauss potential 

As a next example we explore the applicability of different approximations introduced in subsection 3.2. 
For the analysis we use the prescribed potential of Gauss-form: 

Vc<(r) = -2exp(-5r 2 ) (25) 

where both distance r and energy are measured in atomic units (au). 

Potential (l2"5j) provides a set of input phase shifts 5°" s listed in table Q] at scattering energy of E = 18 
au (k — 6 au). Using this input set one calculates the CT inverse potential VctW = E qcT(kr) which is 
depicted in figure [5] Note that Vct is the same as Vq, within the width of line, therefore plot of Vg is 
omitted in figure [2] 
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Figure 1: Effective 87 Rb+ 87 Rb inverse potential Vct (in mK) as a function of the radial distance r (in 
nm) at cm. energy E — 303 /xK. 



The different approximations Va, Vl, and Vr can thus be compared to Vct which may be taken 
to be exact. In figure [2] we see that, for this particular example, the potential VA.( r ) = E q^(kr) is a 
good approximation to the original one in view of its initial behaviour (depth) and asymptotical property 
(range) . Recall that approximation Va is obtained by simply adding inversion approximations V c and V a 
derived by inverting the separate sets S c and S of phase shifts belonging, respectively, to the sets of even 
and odd angular momenta Vs. Approximation Vr is obtained by unifying the calculated sets T c and T Q 
of shifted angular momenta, L's, listed in table Q] under heading (fT7|) or (fT8|) . Finally, the approximation 
Vl has been obtained by simply using the one-term approximate values L a , equation (|24|) . which are also 
listed in tableQ]under heading (f2~4"|). It is interesting to see in figure[5]that this (totally analytical) version, 
the approximation Vl provides a somewhat better result than that of the more involved approximation 
V T . 

Table 1: Original phase shifts 5° r16 produced by the Gauss potential (j25|) at scattering energy E = 18 
au. Shifted angular momenta L corresponding to solution of equations dU), (jl7)) . or (|T8|) . and (|24)) . 
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Figure 2: Inverse potentials V(r) obtained from input phase shifts S°" 6 [see table [T] as a function of the 
radial distance r at energy E = 18 (k = 6) au. Curves obtained by the CT, and approximate methods 
are labeled according to the procedures discussed in the text. 

Table 2: Original phase shifts <5° llg and elasticities rj° Yls taken from the phase shift data of [5] for n 
scattering by 12 C nucleus at the energy E^ ab = 10 MeV. Real (left column) and imaginary (right column) 
part of the shifted angular momenta L listed under headings CT, A, and L, respectively, correspond to 
solutions of equations ©, dTTJ), or (0, and flUD - 
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§ ong 


orig 
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CT 


CT 


A 


A 


L 


L 





0.827 


0.580 


-0.581 


-0.085 
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0.723 
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-0.122 
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-0.365 
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0.021 
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-0.013 



4.3 n + 12 C scattering at E = 10 MeV 

Our last example is taken from nuclear physics. Chen and Thornow [8 have derived 88 sets of complex- 
valued phase shifts from a comprehensive analysis of n scattering by 12 C target nucleus in the scattering 
energy region 7 MeV < E\f < 24 MeV. Because of the spin-orbit coupling, each partial wave provides 
two phase shifts, 8^ and 5~[ . In case of weak spin-orbit coupling the combined phase shift <5/ = [{I + 
1)5 ; + + IS^]/(21 + 1) are characteristic of the underlying central potential [5]. 

One set of such combined phase shifts, valid to the neutron scattering by 12 C at the energy of £7^ = 10 
MeV is listed in table [5J Here, <5° rlg denotes the real part of the combined phase shifts, Re<5i, and ??° rlg 
stands for the elasticity, r]° ng — \ exp(2i<5;)|. The results of the various inversion procedures are shown in 
figure [3] 

In order to draw conclusion about the applicability of the different approximate inverse procedures 
A, T, and L, one can compare the corresponding potential curves to the one obtained by the CT method 
which is assumed the best. (There is no model potential in this case.) Indeed, if one calculates back 
the phase shifts and elasticities that the different inverse potentials Vct, Va, Vt, and Vl are providing, 
we see that the reproduction is the best for the CT potential. All these facts can be studied in table [3] 
where we list the differences A; and between the original data <5° ng and r)° Ils and the re-calculated 
ones provided by the inverse potentials of the respective methods. The potential curves Vct, Va, Vt, and 
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Figure 3: Inverse potentials V(r) obtained from input phase shifts <5°" s and elasticities r]°" s [see table 
H] as a function of the radial distance r at energy E^ h = 10 McV (E°- m = 9.23 MeV, k = 0.638 fm" 1 ). 
Curves obtained by the CT, and approximate methods are labeled according to the procedures discussed 
in the text. The imaginary part of the potentials concentrates around the abscissa, meaning a weak 
absorption. 

Table 3: Differences A; and S/ between the original data and the calculated ones provided by the inverse 
potentials (shown in figure [3]) obtained by the different inverse methods CT, A, T, and L. 
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VL are drawn in figure [3] and we see that the curves obtained by the methods A, T, and L approximate 
well the inverse potential Vqt considered to be the reference. Because the numerical reproduction is also 
surprisingly good we may conclude that the proposed approximations to the CT method can be used for 
a global orientation about the nature of the underlying interaction, if input phase shifts are known from 
another source and the solution of the generic equation (|9]) is not possible. 

5 Conclusion 

By observing the simple analytical property f|13[) of the asymptotic expansion functions A\ (x) of the 
transformation kernel K(x, y) involved in the Cox-Thompson (CT) inverse scattering method at fixed 
energy, simplified solutions of the CT method have been derived [see equations (IT71) and (|T8"|) ] applicable 
to certain special scattering of identical bosons (fermions) when only even (odd) partial waves are arising 
(and thus measurable). 

The new formulas (TT7)) and (fT8")l obtained are easier to be solved by usual nonlinear solvers (such 
as that based, e.g., on the Newton- Raphson procedure [5]). This is because the new formulas do not 
involve the inverse of a matrix containing the unknowns themselves; they 'only' contain product and 
tangent (cotangent) operations. Therefore the range of applicability of the new formulas is wider than 
that of the generic equation ©. However, identical bosonic (fermionic) collision experiments are rare 
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in practice. Therefore we have developed several approximations too, based on the new formulas (fT7|) 
and (|18p . classified according to the level in which they have been introduced. By taking examples from 
atomic collision experiment, potential scattering and nuclear physics, we have demonstrated the wide 
applicability of the new equations which make the solution of the CT inverse scattering method at fixed 
energy easier. 
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